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1, INTRODUCTION 

In this era, the major concern in power system are to discover the power quality due to the impact of 
PQ disturbances (PQDs). The monitoring and analysis of PQD are generally performed with specialized 
measuring equipment based on the standards. Based on the research [1, 3, 5], most of the PQDs where 
generalized in wavelet transform before extracting to the next stage. During the generalized, the minimum 
frequency is fixed in separation threshold in prior and generates the substantial redundant information. 

Moreover, most of the research paper [2, 3, 6], for classify the PQDs it use support vector machine 
(SVM) after the extraction process of the data because the SVM are efficient in large classification problems 
and better in generalization properties compared to conventional classifiers but in others [1, 4, 5], where 
using the classifiers of data which are, decision tree, perception tree and also, random forest that consist a 
group of tree that used. It known as considering the standards and practical signal characteristics. Other than 
that, there are method proposed which contains the feature extraction (FE) and decision-making [7]. In FE, 
histogram are utilize and combine with a Discrete Wavelet Transform (DWT) that can detect the PQ while in 
the stage of decision, Extreme Learning Machine (ELM) classify the dataset of PQ which resultant of the 
performance detected is high [7]. 

The Adaptive Network-based Fuzzy Inference System (ANFIS) introduced by Jang [13, 14] is a 
well-known neural fuzzy controller with fuzzy inference capability had been implemented in various works 
[8-12]. ANFIS is based on fuzzy logic modeling and use ANN as the learning algorithm. 

The aim of this work is to highlight the application of the ANFIS to perform fault disturbances 
classification using three set feature vectors extracted from the S-transform and to find the right type and 
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number of Membership Functions (MF’s) for used in ANFIS. This is to provide a fast action to supply power 
to the consumer if there is some disturbances occur in the power system. 

The remaining parts are organized as follows. Detailed methodology is discussed in Section 2 
Section 3 describes the results obtained. Finally, conclusions are given in Section 4. 


2. RESEARCH METHOD 

The flow chart of the proposed method is shown in Figure |. In this study, three feature vectors 
extracted from the S-transform introduced in [15] were applied as inputs and classification stage with ANFIS 
implementation will be carried out using MATLAB software. 


Inputs: 
Extracted data using 
S-transform method [12] 


Classify and Analyze all fault 
disturbances using ANFIS 





Figure 1. The flow of the proposed method 


2.1. S-transform Feature Extraction 

S-transform is one of the feature extraction methods based on time-frequency analysis. It extracts 
the pre-processing data into the most salient features that represent the power quality phenomenon [16]. 
The original formula for s-transform is shown in (1). 


(r—-t)* Ff? 


st f=- [xe 2 OF ae. (1) 


where: f is frequency, t and T is time 

The output from S-transform is am Xn matrix known as S matrix with column is referred to time 
and row is referred to frequency. Every element from this matrix is a complex value. Each row indicates the 
magnitude of S-transform with all the frequency at the same time and each column indicates the magnitude 
of S-transform from same frequency. 

According to [15], the first step in implementing the feature extraction 1s to map the distorted data 
signal to S-transform domain. For non-stationary signal interference and noise, using the analysis of S- 
transform, feature vectors extracted from signal interference in the form of time-frequency representation of 
the curves. 

A statistical technique is used to the amplitude of contour matrix of S-transform by using maximum 
amplitude and frequency amplitude plot. The resulting features are formulated as follows. 

Amplitude factor, Fl as the first feature with the given by (2), 


fF, = 1 + std] + std2 - normI1- norm2 (2) 


Where: 

std1: Maximum value of standard deviation distortion signal 
std2: Minimum value of standard deviation distortion signal 
norm1!: Maximum value of normal signal 

norm2: Minimum value of normal signal 

Second feature is given by the equation of, 


F2 = stdl —norml (3) 
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And the third feature is given by equation, 
F3 = mean (mean (abs (ds)7)) (4) 


Where: 
ds: absolute value of S-transform distortion signal 

In this study, based on these formulas, three feature vectors extracted were used as inputs in which 
parameter Fl was assigned as input one, input two was represented by parameter F2, and parameter F3 
represented input three. Table 1 depicts the sample of three inputs. 


Table 1. Sample of The Input Data Extracted From S-Transform Method 


eature Vectors 










Types of Fl F2 F3 


Faults 






Sag 1.0573 0.0302 0.0009 
Swell 1.0814 0.045 0.0018 
Impulsive 1.0936 0.059 0.0013 
Interruption 1.0449 0.0386 0.0009 
Notch 1.0703 0.0355 0.0013 
Pure 1.0742 0.0371 0.0014 


2.2. Adaptive Neuro-Fuzzy Inference System 

ANEFIS technique has been considered in this work because it had been stated as one of the simplest 
and has a well-defined mathematical model technique to determine the classification of fault [17-19]. It isa 
hybrid combination of Adaptive Neural Network (ANN) and Fuzzy Inference System (FIS) [12]. FIS use a 
membership mapping modeling of input to determine its output by a membership function parameters. 
This membership function parameter can be adjusted to get the targeted values. ANN function as a learning 
algorithm to the system which gives fast convergence and more accurate in ANFIS target. ANN is adapted so 
that it can computes the membership function parameters that best allow the associated FIS to track the given 
input/output data since the FIS only applied to only modeling systems whose rule structure is essentially 
predetermined by the user's interpretation [12]. 
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Figure 2. Architecture of ANFIS 


addition, ANFIS is applied in this work because it give a high precision, easy, and effective technique 
[12, 13]. 

The architecture of ANFIS is depicted in Figure 2 [20]. It consists of five layers with two inputs and 
one output based on the Takagi-Sugeno model. 
Layer 1: Fuzzification 

This layer is located at the first nodes at the receiving of input data or signal. It adjusted the input 
signal and performed a Membership Function (MF) in (5). The MF parameters are adjusted according to the 
user interpretation. 


Qi = ka(x) (5) 
where ajis membership function of k,(x) and a is linguistic label associated with the label. 
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Layer 2: IF-THEN Stage Rule 
In this stage, the signal that has through a fuzzification from the first layer will represent as the 
firing strength of the IF-THEN rule. The output is multiple of inputs and it performs the fuzzy AND rules. 


W; = Rai(x)*Roi(vy) i= 1,23... (6) 


Layer 3: Normalization 

The input for normalization layer is the output from the IF-THEN rule layer. In this layer, the 
normalized firing strength will be calculated. Every gained signal is divided by the total of gained signal 
using below equation. 


Ri 
— on 
par R 





(7) 


i 


Layer 4: Defuzzification 
In this layer, parameters of output are adjusted because of each node is an adaptive node. 
Normalized signal is gained by the linear function of MF to the output signal. 


fi=Ni(px+qyt+n) i=1,2,3... (8) 


pix, qiy and r; are the MF parameters for the linear signal. 
Layer 5: Neuron Addition 

This is the last layer for ANFIS. In this layer it will calculate the entire defuzzification signal since 
it has only one node, with the summation of the input signal. 


OT = din fl (9) 


3. RESULTS AND ANALYSIS 

Hybrid learning algorithm is employed to determine the optimum values of the FIS parameters of 
the Sugeno-type model. Hybrid-learning is the combination of least-squares method and back-propagation. 
In the ANFIS, the training data were used to train (or tune) a fuzzy model. In this work, a collection of 264 
samples divided into two part; 70% for training and 30% for testing data. 

Before training process, several parameters need to be set up first [21]. The numbers of Membership 
Functions (MF’s) are assigned to each input of the ANFIS and were set to two and three MF’s. The number 
of rules produces for two numbers of MF’s is 8-rules and for three numbers of MF’s, rules generated 1s 27- 
rules respectively. To select the best model fit based on the type of MF’s, two membership functions has 
been selected which are: the Gaussian function (gaussmf) and generalized-bell function (gbellmf). The step 
size 1s the length of each gradient transition in the parameter space. In this project, the parameter of step size 
is set as follows: the initial step size value was set to 0.1; the step size-decreasing rate is set equal to 0.9, 
and the step increase rate was 1.1. Epochs or number of iteration is set to 100. Meanwhile, the hybrid 
learning method was employed, and the linear output type is selected. 

Feature saliency provides a mean for choosing the best features for classification [22]. Since the 
feature vectors used as input may have different levels of relevancy, the changes that happen after training 
the MF’s are examined. Figure 3 shows the plots of the MF’s before and after training. Referring to Figure 3, 
changes of the MF were noted for the case of two and three MF’s using the “gbellmf’ and “gaussmf’. 
In overall, for all the cases, most of MF’s for feature 3 were changed after training process. However, for the 
case of using three MF’s with “gbellmf’, MF’s for feature 2 and feature 3 were changed, whereas no change 
was noted in other plots. These results imply that by implementing different numbers of MF’s and choosing 
the right type of MF’s, the classification results can be improved and in the case of the “gaussmf” it gave the 
best results when three MF’s are in used. 

In addition, the Root Mean Square Error (RMSE) for both training and testing were compared and 
tabulated in Table 2. RMSE error is the training and testing error for the analysis of the parameter. The error 
that produce from this analysis is taken to determine which parameter has the best (the least) error for the 
classification of the fault disturbances. It is obvious that by implementing different numbers and types of 
membership functions, the value of each error rate also changed. According to the Table 2, during training 
process, the ANFIS model with two MFs that used “gaussmf’ had better training learning result compared to 
ANFIS with “gbellmf’. The smallest training error recorded was 0.2956. 
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On the other hand, in term of testing process, in the case of the “gaussmf” it gave the best results 
when three MFs are in used. The lowest testing error recorded was 0.0541. Therefore, it can be said that, 
based on the testing result, the ANFIS with three MF’s had a better classification results compared to other 
ANFIS model that used “gbellmf’. 


Table 2. Comparison Of Error Rate For 2 MFs’ AND 3 MF’s With 2 Types Of MF’s 































































































































































































































Parameter of MF’s Number of MF’s 
2 MF’s 3MF’s 
Training error (RMSE) Testing error (RMSE) Training error (RMSE) Testing error (RMSE) 
Gauss 0.2956 1.6235 0.3902 0.0541 
Gbell 0.3237 1.1106 1.3156 0.3383 
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(b) “Gbellmf” 


Figure 3. Plots for 2 MF’s and 3 MF’s (1) Initial MF’s input and (11-1v) Final MF’s for all three inputs 
(F1, F2, F3) using two types of MF’s (a) “Gaussmf’, and (b) “Gbellmf’ 
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4. CONCLUSION 

In this work, it can be concluded that the analysis of fault classification disturbances based ANFIS 
classifier shows that “gaussmf’ with three numbers of MF’s has suitable method for fault disturbances 
classification. For this study, the use of three gaussian MFs (“gaussmf’”’) has obtained the best classification 
results based on the lowest testing error rate. Therefore, this work shows that type of MF’s and numbers of 
MF’s play an important role in the analysis of the ANFIS classification. Furthermore, the ANFIS classifier 
can improve significantly by manipulating the type and number of the MFs. 
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